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' Abstract. We present an algorithm for obtaining the post-Newtonian expansion of 

CN| ■ the asymptotically flat solution to the Einstein-Maxwell equations describing a rigidly 

rotating disc of dust with a constant specific charge. Explicit analytic expressions are 
O ' calculated up to the eighth order. The results are used for a physical discussion of the 

extreme relativistic limiting cases. We identify strong evidence for a transition to an 
j-H ' extreme Kerr-Newman black hole. 
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^ . 1. Introduction 

Since the formulation of classical general relativity by Einstein in 1915, analytic solutions 
have made a large contribution to the understanding of its meanings and proving 
its properties. A crucial role in this process was played by black hole solutions, as 
their description requires the full generality of Einstein's theory. Therefore black hole 
solutions could be expected to exhibit inherent effects of general relativity. Indeed, 
the Schwarzschild solution containing a single mass parameter M was able to explain 
the perihelion precession of Mercury and the bending of light rays around the sun. 
The Kerr solution extended the class of black hole solutions by the parameter of 
angular momentum J, led to the discovery of ergospheres and permitted the rigorous 
investigation of the Lense-Thirring effect. Furthermore, adding a charge parameter Q 
completes the class of stationary black hole solutions by delivering the most general 
Kerr-Newman solution, see [HE]- 

Interestingly, static spherically symmetric ideal fluid configurations have a minimal 
radius, given by | times the corresponding Schwarzschild radius and thus fail to feature 
a continuous connection to black holes [3]. The analytic solution of the problem of a 
rigidly rotating uncharged disc of dust, in contrast, provides an example for a parametric 
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Figure 1: Schematic of the parameter spaces of black holes and charged dust discs linked 
by two limiting cases. 



transition of equilibrium configurations of ordinary matter to an extreme Kerr black hole 
[I] , see also [5] . Rotating fluid rings admit such a black hole limit too [6] [7] . 

On the other hand, the so called electrically counterpoised dust (ECD) 
configurations, a class of static solutions of dust with a particular constant specific 
charge, possess a parametric transition to the extreme Reissner-Nordstrom black hole, 
see [8]. Inspired by these opposite limiting cases, we conjecture that the dust disc 
with a constant specific charge varying from zero to the ECD value has a parametric 
transition to the extreme Kerr-Newman black hole. The parameter spaces of black 
holes and charged dust discs and their connections are illustrated in figure [U In order 
to investigate this assumption and to support the search for the analytic solution of the 
charged disc, a post-Newtonian expansion is performed. We note that in the uncharged 
case, the analytic disc solution [I] was indeed preceded by the partly numerical but 
highly accurate Bardeen- Wagoner expansion [9], which could later be recovered from 
the analytic solution [10] . The charged disc is expected to show novel effects due to the 
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interplay between electromagnetic and gravitational fields. 

2. Formulation of the disc problem 

2.1. Model of matter 

To facilitate analytic calculations, we employed the following simplifying assumptions: 
Firstly, we use the model of dust, i.e. pressureless matter described by the dust 
part of the energy momentum tensor in (j4j). Moreover, we assume a rigid rotation with 
a constant angular velocity Q. The dust we consider here contains the constant specific 
charge e G [—1,1] so that the charge density p e \ and the baryonic mass density ji are 
connected via 

Pel = ep (1) 

(we use Gauss units and c = 1, G = 1). 

Furthermore, we assume axial symmetry and stationarity described by the Killing 
vectors rj and as well as reflectional symmetry w.r.t. an "equatorial plane". 

2.2. Basic field equations 

We start with Einstein's field equations 

Rab ~ T^RQab = 8nT ab (2) 
and the covariant Maxwell equations 

F [ab;c] =0, F ab ;b = 4nf, (3) 

wherein the energy momentum tensor T ab is simply the sum of an electromagnetic part 
and a dust part: 

-,(em) rp(dust) _ ^ ( t? r? c 1 77 r?cd 



r ab = rr + rr v = ^ [F ac F b c - - 9ab F cd F^j + W6 . (4) 

For the purely convective four-current density holds f = p c \u a = epu a . The metric can 
globally be written in terms of Weyl-Lewis-Papapetrou coordinates as: 

ds 2 = r 1 [h {do 2 + dC 2 ) + o 2 d<p 2 ] -f(dt + ad<p) 2 . (5) 

Moreover, a four-potential of the form A a = (0, 0, A v , A t ) can be introduced by 
F ab = A b . a — A a . b = A bta — A atb , which automatically fulfils the homogeneous Maxwell 
equations. The coordinates correspond to the Killing vectors via d t = £ and = 77, so 
that all the five free functions /, h, a, A t and A^ depend only on the coordinates g and 
(. In addition, fi is made up from a surface mass density by 

H = o-( 0) {g)5{(). (6) 

Note that two other definitions of surface mass density will later be used: a = o~(p)h/ f, 
directly linked with the boundary conditions, and the coordinate independent proper 
surface mass density a p = o^pfjh. 
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The coupled nonlinear system of differential equations that emerge can be 
reformulated as a boundary value problem for the equations valid in the (electro-)vacuum 
outside the disc. 

2.3. Ernst equations 

In the absence of matter, the Einstein-Maxwell equations can be reduced to the Ernst 
equations. They are expressed in two additional potentials /3 and b. Using the 
abbreviations A t = —a and = A, these potentials are defined by 

P,q — ~ ( aoc ,C + Ac) > 0,C = ( aa ,e + Ae)i ( 7 ) 

f? £? 

and 

f 2 f 2 
b, e = a>£ + 2 (/3a jB — otfi, e ) , = — a tB + 2 {fia ^ — «Ac) • (8) 

Q Q 

Thereby the complex Ernst potentials £ and $ can be introduced: 

$ = a + i/3, £ = (/-$$)+ ib. (9) 

In these terms the Ernst equations read [11] 

(m£ + $$) A£ = (V£ + 2$V$) • V£, (10) 
(9?£ + $$) A$ = (V£ + 2<JV$) • V$, (11) 

where the behaviour of V and A is analogous to that in Euclidean 3-space, with 
cylindrical coordinates (g, (, tp). The metric function h is eliminated during this process 
of transforming the field equations to the Ernst equations and can be determined by 
integration afterwards. For the purpose of the post-Newtonian expansion, it is helpful 
to express the Ernst equations in terms of the real functions /, a, f3 and b, which leads 
to 

fAf = fA (a 2 ) + fA (/3 2 ) + V/ (V/ - 2a Va - 2/3V/3) 

-V/3 (V/3 + 2aV/3 - 2/3Va) , (12) 

fAb = 2V/V6 + Aaf3 ((Va) 2 - (V/3) 2 ) - 4 (a 2 - f3 2 ) VaV/3 

+2 V/ (aV/3 - /3Va) - 2V6 (aV« + /3V/3) , (13) 

/Aq = V/Va - V/3 (V6 + 2aV/3 - 2/3Va) , (14) 

fA/3 = V/V/3 + Va (V6 + 2aV/3 - 2/3Va) . (15) 

Boundary conditions on the disc 

To describe the influence of the matter, we have to return to the full Einstein-Maxwell 
equations. We evaluate them in the corotating frame with tp' = <p — Qt, where the metric 
retains its form: 



ds 2 = r 1 



ti (&q 2 + dC 2 ) + Q 2 dp' 2 - /' (dt + a'dtp') 2 (16) 
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with the dashed functions related to the original ones by 
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1 + Qa)f = (l-Qa')f and J = J- 
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Integration over a small flat cylinder around a mass element of the disc (very 
analogous to the well known treatment of surface charge densities in electrostatics) 
delivers matching conditions between the areas above and beneath the disc. These 
matching conditions are transformed to the following four boundary conditions by the 
use of reflectional symmetry: 



(3' = 0, 



v = o, 



a 



« if 



r 



Back in the nonrotating frame they take the form: 

1 + Qa)f 2 
tig 2 



ttg 



[l + Qa)f 



— a ) = A t £, 
= a ,0 



(24) 

(25) 
(26) 
(27) 
(28) 



The treatment in the nonrotating frame is justified by the simple boundary conditions 
at spatial infinity, which are given by g a b — > r]^ and A a — > in that system. 

Overall we have the 4 differential equations (|T2|) - (|T5|) and the 4 boundary conditions 
for the 6 expansion functions /, a, (3, b, A and a. However, due to the 
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Figure 2: Parameter space of the charged disc. 



definitions (J7|) and (jHJ), only 4 of these functions are independent and therefore the disc 
problem is well posed. 

3. Ansatz for the expansion 

3.1. Parameter space of the charged disc 

We introduce the relativity parameter 7, which was already successfully used for the 
post-Newtonian expansion of the uncharged disc [9]: 

1 = 1- VK with / c = /( e = 0,C = 0). (29) 

The parameter 7 is closely related to the redshift Z c of a photon emitted at the center 
of the disc and measured at infinity: 

1 = z c (\ + z c y l . (30) 

Therefore, it is not surprising that the transition to a black hole happens at 7 — > 1 for 
the uncharged and the ECD case. The parameter space of the charged disc as shown 
in figured] (w. l.o. g. restricted to positive charges) is now strongly conjectured to be the 
area 7 G [0, 1], e G [0, 1] shown in figure [2j As a third parameter, the disc radius go is 
added. It naturally comes into play by introducing elliptic coordinates through 

( = Q £rj and g = g ^(l - r] 2 )(l + £ 2 ) 

and carries the whole information about the scale of the disc. By the usage of elliptic 
coordinates and by multiplying powers of £0 to unit carrying quantities, the whole disc 
problem can be made dimensionless. Those normalized quantities, labeled with a star, 
are for example 

a* = a/g , A* = A/g , Q* = g Q and a* = g a. (31) 
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3.2. Formulation of power series 

Inspired by the Newtonian disc of charged dust (in the Newtonian limit we have 
7 = — Un(q = 0, £ = 0) = Q* 2 / (1 — e 2 ) with the Newtonian potential C/n), the squared 
angular velocity is expanded as a power series in 7: 

n* 2 = (1 - e 2 ) 7 + ( 7 2 ) . (32) 

The angular velocity f2* is therefore described by an odd series in 7^: 



DC 



n=l 

All the previously discussed functions are now either symmetric or antisymmetric 
with respect to a change of the sense of rotation Q* — > —Q*. As such, they can be 
expressed either in even or odd power series of 72. The lowest occuring power is 
determined by the Newtonian limit. On the whole we get 



00 00 

3 



n=l n=l n=l (QA\ 

oo 00 00 V / 

n=l n=l n=l 

The remaining coefficients are still functions of £ and rj. These expansions are inserted 
in the 12 equations ( P]) -(fl5]h (125l)-(p5I) plus (0) and ©. These 12 equations, which we 
call expansion equations, are now evaluated by equating coefficients of powers of 7 in 
order to obtain analytic expressions. 

4. Algorithmic solution 

4-1- Structural examination 

The expansion equations in k-th order read: 



(1 + £ 2 )/\ e = -A k Jg + F Dlk 


(35 a) 


(1 - r] 2 )(3 k , v = +A kti /g + F D2k 


(356) 


(1 + £ 2 )&fc,£ = +a k)rt /g + F D3k 


(35 c) 


(1 - ri 2 )b k , n = -a k ^/g + F D4k 


(35 d) 


Vl - e 2 (l - ri 2 ) a k,i = +A k ^/ g + F Blk 


(36 a) 


2a/1 -e 2 (l -ri 2 )fk,( = -a k ,t/g + F B2k 


(366) 


e fk,( = — 2a fc ,5 + -FB3k 


(36 c) 


Wl - e 2 VQo^k + /*,,, = -2ea M + F B4 k 


(36 d) 


A/ fc = F Elk 


(37a) 


A6fe = F E 2k 


(376) 


Aa fc = F E3 k 


(37c) 


A/3fc = F E 4k 


(37 d) 
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For clarity, only coefficient functions of k-th order have been written down, and terms 
of lower orders in each equation are subsumed in the Symbol F. These functions F all 
vanish for k = 1 so that the first order equations reduce to coupled Laplace equations 
with boundary conditions. The gravitational and the electric potential of the Newtonian 
disc are reproduced through fi = 2U^ and a± = —U c \ (see ( 1A.2I) ). For k > 1, none of 
the functions F vanish and we get coupled Poisson equations. The requirement for 
asymptotical flatness implies that all the 6 coefficient functions f k , cik, At, b k , A k and 
a k have to vanish at spatial infinity. 

4-2. Solving the boundary value problem 

A very helpful observation comes from the fact that all expansion functions are 
polynomials in rj. Therefore the functions F (besides a factor of (£ 2 + r/ 2 ) -1 in the 
-^Eik) are also polynomials in r\ and the Poisson inhomogeneity can thus be fragmented 
to Legendre polynomials. Then the Poisson equation with an inhomogeneity consisting 
of a single Legendre polynomial, 

^ = { [d + {"Me] f + [(1 - ,„} = '-f^, (38) 

is transformed by the ansatz ip = A(!;)P n (rj) into the ordinary differential equation 

[(l + £ 2 )A^-n(n + l)A = /(£). (39) 
Two independent (real) solutions to the corresponding homogeneous equation are 

A x = i n P„(i£), A 2 = i 1 ~ n g„(iO (40) 

where Q n are the Legendre functions of the second kind. Now the inhomogeneous 
solution can be composed as 

A( « = - M * II TT^W Ax + Mi) II ^W Ax ' w = A '° Al - A[A - (41) 

With adequately chosen constants ci and C2, the boundary conditions at C, — > oo and 
£ — > can be fulfilled respectively. 

4-3. Completion of the expansion equations and algorithmic solution 

In addition to the 12 expansion equations f)35ap - (137d|) . two parameter relations are 
necessary to determine an integration constant and the scalar component for each 
order. First, from ( 129]) . we immediately obtain at g = 0, £ = 

Ac = -2, /ac = 1 and f kc = for fc > 2. (42) 

Secondly, the following relation to the surface mass density can be derived analogously 
to the derivation of the boundary conditions (F___ still collects lower order terms): 

- [f k £ + -FREGlk] = 87T^ O"fe + -pREG2k (43) 
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Figure 3: Plot of Q*, the angular velocity normalized by the disc radius. 

In order to have a finite value of o~k at the rim of the disc, rj — 0, the expression 
[fk,£ + -pREGik] has to contain a global factor r\. The integration constants of the other 
unknown functions vanish due to reflectional symmetry. 

The system of expansion equations can be dealt with by setting up the linear 
combination O k = ~e/& + ct k . This provides a quantity with a given Poisson 
inhomogeneity and a given boundary condition: 

AO fc = |F Elk + F E3k , O u = ^F B3k . (44) 

Once O k is determined according to 4.2, f k can be determined by (136 dj) and a k is 
immediately obtained. Hereupon (3k and bk can be computed successively out of Poisson 
boundary value problems and A k and a k eventually by integration. 

5. Results 

With the methods described above and making use of Mathematica by Wolfram 
Research, coefficient functions of the 6 expansion functions could be calculated up to 
the eighth order. The first important observation is that these coefficient functions are 
not only polynomial in r\ but also in £ and arccot £. The first 2 orders are listed in the 
appendix, higher orders are available as data files. Noteworthy are the global prefactors 
in the expansion functions, which could be found in all of their coefficients: 

Q* k oc Vl - e 2 , a k oc e, (3 k oc ey/l- e 2 r], b k oc \/l - e 2 r], 

A* k (xeVT^(l-ri 2 ) and a* k oc v / T 3 ^(l - r, 2 ). 

Within the expansion, the angular velocity Q* is also calculated. It is depicted in figure [3] 
up to the eighth order. Here the parameter space (figure [2]) is shown in the 7-e plane. 
Q* vanishes at the transition to the static configurations for e — > 1 and in the extreme 
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relativistic limit (7 — > 1). The latter is caused by the vanishing disc radius qq. This 
shows that £0 is no longer a well chosen parameter for 7 — >■ 1. 

The further evaluation of the results is done by discussing global quantities of the 
disc. The gravitational mass M, the charge Q and the angular momentum J can be 
conveniently deduced from the far field behaviour of the expansion functions: 

f= i-^|- + o(r 2 ), (46) 



Q* 

a = T 



O (r 2 ) , (47) 



2\ 



One compelling feature of the disc with constant specific charge, is that even the baryonic 
mass Mo = e~ 1 Q can be derived from the asymptotical behaviour. Based on this, the 
relative binding energy = (M — M)M _1 can be directly calculated (even in the 

limit e — > 0). It is depicted in figure H] up to the seventh significant order. The first 
three orders read 

<" = i (1 - O 7 + ± - 9^ + 7) * + 3^ [179200 ( 3e * - ^ + 64) 

-7T 2 (58618e 4 - 1042821e 2 + 1275648)] 7 3 + O (7*) , (49) 

all seven orders are given in the appendix. At this point, the accuracy of the expansion 
can be tested by comparison with the known analytic solution of the uncharged disc. 
The potentially most critical point is e — >• 0, 7 — >• 1, where the transition to the 
extreme Kerr black hole occurs. Here the numerically evaluated expansion gives 
E^ el) (e = 0,7 = 1) = 0.3614, which differs from the analytic value 0.373283588... 
[5] by about 3%. In the uncharged limit (e — > 0), the expressions for the remaining 
nontrivial functions / and a agree with an expansion of the analytic solution of the 
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uncharged disc |TU]. In addition, the results for the global quantities for e — > comply 
with those of the uncharged disc given in [5]. A non-trivial test for all e consists in 
checking the exact relation [12] 

M = 2QJ + (1 — 7 + ea c )M , a c = a{g = 0, ( = 0). (50) 

This relation is satisfied by our results for all eight available orders in 7. 

Figure [5] shows the coordinate independent surface mass density a p normalized by 
the angular velocity Q. In contrast to the normalization by £0 i n figure [31 this gives a 
plot regular at 7 — > 1. The already known shift of the maximum of a p /fi out of the 
disc centre appears to weaken with growing charge. For e — > 1 the angular velocity Q 
vanishes and is hence no more an appropriate scaling parameter. 




Figure 5: Surface mass density normalized by angular velocity, cr p /fi, for specific charges 
of (a) e = 0, (b) e = 0.3, (c) e = 0.7, (d) e = 0.99. 

Another interesting quantity is the disc radius go, normalized by the baryonic mass 
Mo which is depicted in figure [6j The small values at 7 — > 1 and the weak dependence 
on e indicate a transition to a black hole as in the limiting cases e = and e = 1. Note 
that one has to distinguish the external perspective (finite g/M, C/M) from the internal 
perspective (finite g/go, (/go) in the limit 7— > 1, see the discussion in [12] . 

6. Conclusions 

The vanishing of the quantity go/M for 7 — > 1 can be interpreted as the contraction 
of a disc with given baryonic mass to the origin of the {g/M, C/M) -coordinate system 
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Figure 6: Disc radius go normalized by the baryonic mass M . 



in the extreme relativistic limit. This strongly indicates the occurrence of a singularity 
in that limiting case. To check whether this singularity is indeed related to an extreme 
Kerr-Newman black hole, three other, more specific tests are carried out. For the 
extreme Kerr-Newman metric, the following relation holds on the horizon Ti (located 
at q/M = = (/M): 

-l 



a\H) 



M 



Q_ 

M 



(51) 



The corresponding e-dependence of a 1 '(q/M = 0, (/M = 0) computed out of 



(rel) 



via M/Q 



'l-E- 



(relh 



is depicted in figure [7J as the red line. If we transform to 



normalized coordinates by g — > g*, ( — > (keep in mind that we expect go for 
7 — >■ 1), the extreme relativistic limit yields the limiting spacetime from the internal 
perspective. The relation a'(g*) = const on the whole disc is strongly conjectured to be 
necessary and sufficient for obtaining the extreme Kerr-Newman metric as the external 
solution (i.e. for (g 2 + ( 2 )/M 2 > 0), provided g — > as 7 — > 1. The curves of a'(g* = 0) 
and a'(g* = 1) are depicted in figure [7] in blue and green. The similarity of the three 
curves shows that a' is indeed in good approximation independent of g* and that its 
value meets the one given by ( 15T|) . 

Furthermore the extreme Kerr-Newman black hole obeys the 2 parameter relations 

J/M 2 



■ M 2 t" M 4 



and 



(52) 



OM(l + J 2 /M 4 ) 

with Q meaning the angular velocity of the horizon. The first of these quotients, 
calculated out of the disc quantities, is shown in figure [HJ For the generic disc Q\ 
is greater than 1 as expected, whereas in the extreme relativistic limit Qi converges 
to unity for growing expansion order n. The last statement also holds for the second 
quotient Q2 shown in figure [9j Therefore we conclude that we have found strong evidence 
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Figure 7: a' calculated via ()5ip (red) as well as taken from the disc values at g* = 
(blue) and g* = 1 (green) in the limit 7—^1 (expansion order n = 8). 




Post-Newtonian expansion of a disc of charged dust 



14 



1.15 - 



1.10 



1.05 



1.00 




0.0 0.2 0.4 0.6 0.8 1.0 

e 

Figure 9: Plots of Q2 in the limit 7 — > 1 for different expansion orders n. 

for the transition of the charged disc to an extreme Kerr- Newman black hole in the limit 
7 -> 1. 
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Appendix A. Coefficient functions in first order 

In the appendix we use the abbreviation X := arccot£. Global factors are exposed at 
the cost of showing the fragmentation to Legendre polynomials. 
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Appendix B. Coefficient functions in second order 
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h = gg-2 [(-72tt + 157re 2 ) £ + (72 + 36e 2 ) £ 2 + (-108tt + 1537re 2 ) £ 3 

+ [-72tt + 157re 2 + (144 + 72e 2 ) £ + (108tt - 667re 2 ) £ 2 

+ (-144 - 72e 2 ) £ 3 + (108tt - 1537re 2 ) £ 4 ] X + [72 + 36e 2 

+ (-144 - 72e 2 ) £ 2 + (72 + 36e 2 ) £ 4 ] X 2 + rj 1 { (612tt - 6067re 2 ) £ 

+ (-432 - 216e 2 ) £ 2 + (1080tt - 15307re 2 ) £ 3 + [-108tt + 66ne 2 

+ (-288 - 144e 2 ) f + (-972tt + 11167re 2 ) £ 2 + (864 + 432e 2 ) £ 3 

+ (-10807T + 15307re 2 ) £ 4 ] X + [144 + 72e 2 + (288 + 144e 2 ) £ 2 

+ (-432 - 216e 2 ) £ 4 ] X 2 } + V 4 { (-660tt + 9357re 2 ) £ 

+ (648 + 324e 2 ) { 2 + (-1260tt + 17857re 2 ) { 3 + [108tt - 1537re 2 

+ (-432 - 216e 2 ) f + (l080vr - 15307re 2 ) £ 2 + (-1296 - 648e 2 ) £ 3 

+ (1260tt - 17857re 2 ) £ 4 ] X + [72 + 36e 2 + (432 + 216e 2 ) £ 2 

+ (648 + 324e 2 )£ 4 ]X 2 }] , (B.2) 

a 2 = [(81tt - 247re 2 ) £ - 72£ 2 + (27vr - 727re 2 ) £ 3 + [81tt - 247re 2 - 144£ 

+ (-90tt + 487re 2 ) £ 2 + 144£ 3 + (-27tt + 727re 2 ) £ 4 ] X + [-72 + 144£ 2 

-72£ 4 ] X 2 + 7] 4 { (165tt - 440vre 2 ) £ - 648£ 2 + (315tt - 8407re 2 ) £ 3 

+ [-27tt + 727re 2 + 432£ + (-270tt + 7207re 2 ) £ 2 + 1296{ 3 + (-315vr 

+840vre 2 ) £ 4 ] X + [-72 - 432£ 2 - 648£ 4 ] X 2 } + r/ 2 {(-342tt 

+336vre 2 ) f + 432£ 2 + (-270vr + 7207re 2 ) £ 3 + [90tt - 48vre 2 + 288£ 

+ (432tt - 5767re 2 ) £ 2 - 864£ 3 + (270tt - 7207re 2 ) £ 4 ] X 

+ (-144- 288£ 2 + 432£ 4 ] X 2 }] , (B.3) 



Pi 2 

P 2 = X44(hr2 71 [( 9225?r - 55207re 2 ) { 2 - 12960£ 3 + (-7425vr - 7200vre 2 ) £ 4 

+ [(-37357T + 12007re 2 ) f - 8640£ 2 + (-6750tt + 79207re 2 ) £ 3 + 25920£ 4 
+ (7425tt + 7200vre 2 ) £ 5 ] X + [4320{ + 8640£ 3 - 12960£ 5 ] X 2 
+r/ 4 {-21127T - 2048vre 2 - 17280£ + (-24255tt - 235207re 2 ) { 2 - 64800£ 3 
+ (-311857T - 302407re 2 ) £ 4 + [5760 + (7425vr + 72007re 2 ) £ + 77760£ 2 
+ (34650vr + 33600vre 2 ) £ 3 + 129600£ 4 + (31185tt + 302407re 2 ) £ 5 ] X 
+ [-12960£ - 60480£ 3 - 64800£ 5 ] X 2 } + rf 2 {-4320tt + 2240vre 2 + 5760£ 
+ (1950tt + 260007re 2 ) £ 2 + 60480£ 3 + (34650tt + 336007re 2 ) £ 4 + [5760 
+ (6750tt - 79207re 2 ) f - 23040£ 2 + (-13500tt - 372007re 2 ) £ 3 - 120960£ 4 
+ (-346507T - 336007re 2 ) £ 5 ] X 

+ [-8640£ + 17280£ 3 + 60480£ 5 ] X 2 }] , (B.4) 



VT-e 

1447T 2 



2 

-r] [(-3024tt + 27037re 2 ) £ 2 + 3456£ 3 + 38257re 2 £ 4 + [(1008tt - 56l7re 2 ) £ 
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+2304£ 2 + (3024tt - 39787re 2 ) £ 3 - 6912£ 4 - 38257re 2 £ 5 ] X + [-1152f 

-2304£ 3 + 3456£ 5 ] X 2 + rj 2 {l344vr - 10887re 2 - 1536£ + (5040tt 

-13430vre 2 ) £ 2 - 16128£ 3 - 178507re 2 £ 4 + [-1536 + (-3024tt 

+3978vre 2 ) £ + 6144£ 2 + (-5040tt + 193807re 2 ) £ 3 + 32256£ 4 

+ 178507re 2 ^ 5 ] X + [2304£ - 4608£ 3 - 16128£ 5 ] X 2 } + rj A {l088vre 2 

+4608£ + 124957re 2 £ 2 + 17280£ 3 + 160657re 2 £ 4 + [-1536 - 38257re 2 { 

-20736£ 2 - 178507re 2 £ 3 - 34560£ 4 - 160657re 2 £ 5 ] X 

+ [3456{ + 16128£ 3 + 17280£ 5 ] X 2 }] , (B.5) 

A* 2 = 57 ~J (1 - V 2 ) [-1536 + (423tt - 727re 2 ) £ + 288£ 2 + (840vr - 6327re 2 ) £ 3 
-864£ 4 + (-495tt - 4807re 2 ) £ 5 + [423vr - 72vre 2 + 576£ + (-747tt 
+240vre 2 ) £ 2 + (-675tt + 7927re 2 ) £ 4 + 1728£ 5 + (495vr + 480vre 2 ) £ 6 ] X 
+ [288 + 864£ 2 - 288£ 4 - 864£ 6 ] X 2 + rf {-1536 + (-3729vr 
-36167re 2 ) £ - 19296£ 2 + (-13860tt - 13440vre 2 ) ( A - 21600£ 4 
+ (-103957T - 100807re 2 ) £ 5 + [495tt + 480vre 2 + 12096£ + (7425tt 
+7200vre 2 ) £ 2 + 52992£ 3 + (17325tt + 168007re 2 ) £ 4 + 43200£ 5 + (10395vr 
+ 10080vre 2 ) £ 6 ] X + [-864 - 12960£ 2 - 33696£ 4 - 21600£ 6 ] X 2 } 
+rj 2 {-1536 + (-4014vr + 23767re 2 ) £ + 4032£ 2 + (1740tt + 89207re 2 ) £ 3 
+ 12096£ 4 + (6930tt + 67207re 2 ) £ 5 + [1170tt - 3127re 2 + 3456£ + (4050vr 
-4752vre 2 ) £ 2 - 16128£ 3 + (-4050tt - 111607re 2 ) £ 4 - 24192£ 5 
+ (-69307T - 67207re 2 ) £ 6 ] X 

+ [-576 - 576£ 2 + 12096£ 4 + 12096£ 6 ] X 2 }] , (B.6) 

a* 2 = v (l - V 2 ) [-2048 + 512e 2 + (504tt - 1537re 2 ) £ + (1152 - 288e 2 ) £ 2 

+ (1512tt - 1564vre 2 ) £ 3 + (1152 - 288e 2 ) £ 4 - 12757re 2 £ 5 + [504vr 
-1537re 2 + (2304 - 576e 2 ) £ + (-1008tt + 56l7re 2 ) £ 2 + (— 1512tt 
+1989vre 2 ) £ 4 + (-2304 + 576e 2 ) £ 5 + 12757re 2 £ 6 ] X + [1152 - 288e 2 
+ (-1152 + 288e 2 ) £ 2 + (-1152 + 288e 2 ) £ 4 + (1152 - 288e 2 ) £ 6 ] X 2 
+r] 4 {-2048 + 512e 2 - 96057re 2 £ + (4224 - 1056e 2 ) £ 2 - 357007re 2 £ 3 
+ (5760 - 1440e 2 ) £ 4 - 26775vre 2 £ 5 + [l275vre 2 + (-768 + 192e 2 ) f 
+ 191257re 2 ^ 2 + (-12288 + 3072e 2 ) £ 3 + 44625vre 2 £ 4 + (-11520 
+2880e 2 ) £ 5 + 26775vre 2 { 6 ] X + [1152 - 288e 2 + (3456 - 864e 2 ) £ 2 
+ (8064 - 2016e 2 ) £ 4 + (5760 - 1440e 2 ) £ 6 ] X 2 } + rj 2 {-2048 + 512e 2 
+ (-65527T + 58147re 2 ) £ + (-6912 + 1728e 2 ) £ 2 + (-7560tt 
+231207re 2 ) £ 3 + (-2304 + 576e 2 ) £ 4 + 178507re 2 £ 5 + [1512tt - 714vre 2 
+ (-4608 + 1152e 2 ) £ + (9072vr - 119347re 2 ) £ 2 + (7560vr - 290707re 2 ) £ 4 
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+ (4608 - 1152e 2 ) £ 5 - 178507re 2 £ 6 ] X + [2304 - 576e 2 + (11520 
-2880e 2 ) £ 2 + (6912 - 1728e 2 ) £ 4 + (-2304 + 576e 2 ) £ 6 ) X 2 )] . (B.7) 

Appendix C. Relative binding energy up to the seventh order 

& = (1 - {\, + 4 (7 - 2 £ >) y + [-!79 2 00 (64 - 52^ + 

+7T 2 (1275648 - 1042821e 2 + 58618e 4 )] 7 3 

[-1433600 (-4100096 + 4180767e 2 - 513631e 4 

1192181760000tt 2 L v 

+63360e 6 ) + vr 2 (-619450073088 + 643913271759e 2 - 89452005694e 4 
+8224514048e 6 )] 7 4 

+ r [601896716861440000 (3e 8 - 79e 6 + 556e 4 

3808877661388800000tt 4 L v 

-992e 2 + 512) - 1433600vr 2 (296565800960e 8 - 5672224097024e 6 



+23997689785665e 4 - 40280188824129e 2 + 22326414409728) 
+7T 4 (22758910736859136e 8 - 352585265954279424e 6 
+126049642449134226e 4 + 181545217870667751e 2 
+ 124003773931585536)] 7 5 

+ - [vr 4 (-497795548816521583132672e 10 

157259722042569129984000000vr 4 L v 

+85779761 13717247992135680e 8 - 79165709064202270230097920e 6 

+33090553045688299692921150e 4 + 151318800941422343720072625e 2 

-110706740848531490092351488) - 3064201467658240000 (l7842176e 10 

-386936448e 8 + 3563361915e 6 - 14430944246e 4 + 20418008059e 2 

-9181331456) + 1433600tt 2 (7439589986166374400e 10 

-141373329772495503360e 8 + 1305344205610400847104e 6 

-3328474247978028922343e 4 + 3357975255903479587047e 2 

-1218264268983760846848)] 7 s 

1 

+ 42049186426629786850736209920000000vr 6 ^ 

-23625990069063869264830660608000000 (3e 4 - 52e 2 + 64) (e 4 - 9e 2 + 8) : 
+ 1924145348608007r 2 (l38824602612452556800e 12 
-3668983313237803008000e 10 + 31662317283432188592128e 8 



-128804062549475170301551e 6 + 254210776514571609476318e 4 
-232717498273294565021295e 2 + 79178625735391287705600) 
-204800tt 4 (I4969569145984059227242496000e 12 
-309560598463979551791045935104e 10 
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+1708384495262927033407154421760e s 

-6610208901059681863920531779328e 6 

+10529097859427840958591426694917e 4 

-6453024674363267467191440465925e 2 

+1142467883386445362059993415680) 

+3vr 6 (36492994140024750194896213639168e 12 

-597916823876067055030014970953728e 10 

+1607091649292843235240088719851520e 8 

-15422004663792240785181372461404160e 6 

+19532343512100245042235455926381170e 4 

+5613107997015131353134697962936567e 2 

-10612032886678676548237530771750912)] 7 7 } + O (7 s ) . (C.l) 
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